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ABSTRACT 

We present new observational results that conclude that the nearby radio galaxy 
B2 0722+30 is one of the very few known disc galaxies in the low-redshift Universe 
that host a classical double-lobed radio source. In this paper we use Hi observations, 
deep optical imaging, stellar population synthesis modelling and emission-line diag- 
nostics to study the host galaxy, classify the Active Galactic Nucleus and investigate 
environmental properties under which a radio-loud AGN can occur in this system. 
Typical for spiral galaxies, B2 0722+30 has a regularly rotating gaseous disc through- 
out which star formation occurs. Dust heating by the ongoing star formation is likely 
responsible for the high infra-red luminosity of the system. The optical emission-line 
properties of the central region identify a Low Ionization Nuclear Emission-line Re- 
gion (LINER)-type nucleus with a relatively low [O III] luminosity, in particular when 
compared with the total power of the Fanaroff & Riley type-I radio source that is 
present in this system. This classifies B2 0722+30 as a classical radio galaxy rather 
than a typical Seyfert galaxy. The environment of B2 0722+30 is extremely H I-rich, 
with several nearby interacting galaxies. We argue that a gas-rich interaction involving 
B2 0722+30 is a likely cause for the triggering of the radio- AGN and/or the fact that 
the radio source managed to escape the optical boundaries of the host galaxy. 

Key words: galaxies: active - galaxies: individual: B2 0722+30 - galaxies: interac- 
tions - galaxies: jets - galaxies: stellar content - ISM: kinematics and dynamics 



1 INTRODUCTION 

Radio galaxies are systems that contain an Active Galactic 
Nucleus (AGN) from which powerful radio-continuum jets 
emerge. In many cases, these radio jets reach far beyond the 
optical boundaries of the host galaxy. In the nearby Uni- 
verse, the host galaxies of classical radio sources are gener- 
ally early-type systems. IVeron- Cettv fc Veronl (|200lh argue 
that almost all classical low-z double radio sources are asso- 
ciated with giant elliptical galaxies and that the more discy 
SO hosts are often misidentified ellipticals. In the hierarchi- 
cal model of galaxy formation, early-type g alaxies form the 
end products of merging systems. Indeed, He ckman et al.l 
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l|l986f ) and ISmith fc Heckmanl (|l989f ) showed that a signif- 
icant fraction of the more powerful, edge -brightened radio 
galax ies of Fanaroff & Riley type-II (FR-II: lFanaroff fc Rilevl 
1 1974T ) still show optical signatures of the merger event (such 
as tails, bridges, shells, etc.). Lower power, edge-darkened 
Fanaroff & Riley type-I (FR-I) radio sources are gener- 
ally also hosted by genuine elliptical galaxies, but signs 
of a merge r event in these system s are often much more 
subt le (e.g. Colina fc de Juan 19951) or not apparent at all 
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(e.g. iHeckman et al.lll98rj ; iBaum et al.lll992l : Emo nts et al.l 
l2007| y Recent studies suggest that the difference in host 
galaxy properties between various types of radio sources 
may be more closely related to the excitatio n properties 
of the AGN rather than merely the FR-type (lAllen et al.l 
l200rj : lHardcastle et alj|2007t iBaldi fc Capettill2008j ). but the 
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fact remains that radio galaxies in the nearby Universe are 
typically ear l y-typ e systems. Also for low redshift QSOs, 
ISmith et al.l ( 19861 ) find a trend that the host galaxies of 
radio-loud objects tend to be better fitted by elliptical 
galaxy models and show indications that they are more likely 
to have optical peculiarities than radio-quiet hosts, which 
are better fi tted by disc ga l axy m odels. However, in a more 
recent study iDunlop et al.l (|2003l ) conclude that, for nuclear 
luminosities above the traditional quasar-Seyfert boundary 
of My(nuc) ~ —23, all low-redshift quasars (both radio- loud 
and radio-quiet) are hosted by massive ellipticals. 

In contrast , AGN in spiral galaxies are typically Seyfert 
nuclei l|Sevfertlll943l ), which contain a bright emission-line 
region as a result of an active nucleus that is luminous in 
the optical, ultra-violet and X-ray (though not as luminous 
as typical QSOs/quasars). Seyfert galaxies often also con- 
tain a nuclear radio source, but in contrast to radio galax- 
ies these sources are ge nerally weak and very compact (e.g. 
iMiddelberg et aill2004l . and references therein). Radio- loud 
Seyferts (i.e. with radio powers in the regime of radio galax- 
ies) do exist, but they are rare and their radio compo- 
nent generally does n ot reach beyond the n uclear region of 
the host galaxy (e.g. iKomossa et alj l2006t and references 
therein). 

Although it is very rare for genuine disc galaxies in the 
nearby Universe to host an extended radio source, two clear 
counter-examples exist that have been studied in detail; the 
spiral galaxy 0313-192 in the Abell cluster A4 28 hosts a gi- 
ant ( 3 50 kp c ) double-lobed FR -I radio source l|Ledlow et al.l 
1 19981 . 1200 ll ; (Keel et all [2001 ). while NGC 612 is a typi- 
cal SO galaxy with a large-scale sta r-forming HI disc and 
a powerful FR-I/FR-II ra dio source (|Veron-Cettv fc Veronl 
l200ll ; lEmonts et "all 120081 . see Sect. [O] for more details on 
these two systems) . A third example of a spi ral-hosted radio 
galaxy is the BL Lac obje ct PKS 1413+135 (iMcHardv et all 
Il994l ; lLamer et al.l fl999). but the properties of its radio 
source are difficult to investigate, given that the radio jets 
are pointing roughly in the line-of-sight of the observer. 
These rare cases of disc-dominated radio galaxies provide 
an excellent opportunity for studying host galaxy proper- 
ties and environmental effects that could be important for 
the triggering and/or evolution of the radio source in these 
systems. This could provide a better insight in the origin 
and evolution of radio-loud AGN in general. In addition, 
a detailed knowledge of these nearby disc-dominated radio 
galaxies provides valuable information for comparison stud- 
ies at high redshift, for which recent results suggest that spi- 
ral host galaxies may be muc h more common am ong classical 
radio sources than at low z l|Norris et al.ll2008l ). 

In this paper we study another nearby disc-dominated 
galaxy that hosts an extended radio source, namely 
B2 0722+30. In the literature, B2 0722+30 has occasion- 
ally been mentioned as a disc-dominated radio galaxy, but 
it has not been discussed in detail, mainly because its AGN 
is often confused with that of Seyfert galaxies. In this paper 
we investigate the optical and radio properties of the AGN in 
B2 0722+30 in order to obtain a more reliable classification. 
We also study the morphology, gas content and star forma- 
tion properties of the host galaxy and its environment. This 
allows us to verify whether B2 0722+30 is indeed a rare ex- 
ample of a disc-dominated radio galaxy and investigate the 



conditions under which an extended radio source can occur 
in this system. 

1.1 B2 0722+30 

The radio source B2 0722+30 has a two-sided jet with 
a slightly distor ted structure l|de Ruiter fc Parmal 1 19841 ; 
iFanti et al]|l986h . Being associated to a nearby galaxy (z = 
0.0189), its total radio p ower is P 40 smhz = 1.3x 10 24 W Hz" 1 
(|Giovannini et a l. 20050. Although this is at the low end of 
the radio continuum power for radio-loud AGN, the radio 
structure is 13.6 kpc in extent, reaches outside the optical 
host galaxy and h a s an FR-I morphol ogy l|Parma et aljl986l ; 
IFanti et all 1 19861 ; IFanti et ail Il987f). From H ubble Space 
Telescope (HST) imaging, ICapetti et ail (|2000h claim that 
B2 0722+30 is hosted by a spiral galaxy. Dust lanes are 
clearly visible along the disc and bulge of this system. A 
CO detection traci ng 1.2 x 10 9 My o f Hg g as in B2 0722+30 
w as discovered by iMazzarella et ail (|l993l ) and followed-up 
bv lEvans et~ail (|2005|) using sin gle-dish observations. In a re- 
cent paper l Emonts et alj|2007t l . we presented a preliminary 
analysis of the neutral hydrogen (H I) gas in B2 0722+30, 
which shows that an H I disc follows the optical disc in this 
system. A more detailed H I analysis is given in this paper. 

The outline of the paper is as follows. In Sect. [5] we de- 
scribe the new observations presented in this paper. Section 
[3] investigates the properties of the AGN, the host galaxy 
and the environment of B2 0722+30 (optical morphology, 
gas content, star formation properties and AGN line diag- 
nostics). This is followed in Sect. [4] by a detailed discussion 
about the classification of B2 0722+30 as a disc-dominated 
radio galaxy and the conditions under which an extended 
radio source occurs in this system. In Sect. [5] we give the 
final conclusions. 

For simplicity, in this paper, we will use the name 
B2 0722+30 for both the radio source and the host galaxy. 
We use Ho = 71 km s~ x Mpc -1 throughout this paper, 
which puts B2 0722+30 at a distance of 80 Mpc (1 arcsec = 
0.39 kpc). 



2 OBSERVATIONS 
2.1 Neutral hydrogen 

HI observations of B2 0722+30 were done with the Very 
Large Array in C-configuration on 23 December 2002. The 
effective integration time on the source was 3.1 h. We used 
the 6.25 MHz band with 64 channels in 2 IF mode. For 
the data reduction and visualisation we used MIRIAD and 
KARMA. After flagging and calibration, we fitted a straight 
line to the line-free channels in order to separate the con- 
tinuum from the line data. Even at uniform weighting the 
radio source in B2 0722+30 is only marginally resolved. For 
a detail ed image of t he ra dio continuum of B2 0722+30 we 
refer to IFanti et ail l|l986l . also shown in Fig. Illj l . The line 
data were used to construct a natural weighted data-cube 
with a beam of 19.0 x 17.5 arcsec 2 , velocity resolution of 
43 km s" 1 (after Hanning smoothing) and a noise level of 



1 corrected for Ho = 71 km s 1 Mpc 1 . 
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0.19 mjy beam -1 . From this data set a mask was created 
by smoothing the data spatially by a factor 1.4 and subse- 
quently masking out all the signal below 3er plus discarding 
the signal from the absorption. This mask was used to ex- 
tract a total intensity image of the emission-line gas in the 
data cube by adding the signal in the regions that were not 
masked out (Fig. [2)| . In order to visualise in detail the sig- 
nal in the separate channels (Fig. [3j , another data-cube was 
made with robust +1 weighting and smoothed to a beam- 
size of 26.5 x 24.4 arcsec 2 . The velocity resolution of this 
data cube is 43 km s -1 (after Hanning smoothing) and the 
noise level is 0.23 mjy beam -1 . 



2.2 Optical imaging 

A 3600 sec exposure B-band image of B2 0722+30 and its 
environment was obtained on 13 March 2007 at the 2.4m 
Hiltner telescope of the MDM observatory, located at the 
southwestern ridge of Kitt Peak, Arizona (USA). Observa- 
tions were done under photometric conditions, at an airmass 
of sec z < 1.2 and with a seeing of about 1.1 arcsec. We 
used the Image Reduction and Analysis Facility (IRAF) to 
perform a standard data reduction (bias subtraction, flat- 
fielding, frame alignment and cosmic-ray removal). Likely 
due to minor shutter issues during readout, a gradient was 
present in the background of the 9.5 x 9.5 arcmin Echelle 
CCD image. We were able to remove this effect to a high 
degree by fitting the background in the region surrounding 
B2 0722+30 and subsequently subtracting this background- 
gradient from the image. Using KARMA, we applied a world 
coordinate system to the image by identifying a few dozen of 
the foreground stars in a Sloan Digital Sky Survey (SDSS) 
image of the same region. The newly applied coordinate sys- 
tem agrees with that of the SDSS image to within 1 arcsec, 
good enough for an accurate overlay with our HI data. 



2.3 Optical spectroscopy 

Optical long-slit spectra of B2 0722+30 were taken at the 
William Herschel Telescope (WHT) on 12 January 2004 us- 
ing the ISIS long-slit spectrograph with the 6100 A dichroic, 
the R300B and R316R gratings in the blue and red arm and 
the GG495 blocking filter in the red arm to cut out sec- 
ond order blue light. This resulted in a wavelength coverage 
from about 3500 to 8000 A and A-resolution of about 5 A. 
The slit had a width of 1.3 arcsec and was aligned along the 
major axis of the host galaxy (PA 141°). Observations were 
done at an airmass of about 1.0 and seeing conditions varied 
from 1.2 to 1.6 arcsec during the observations. The total in- 
tegration time was 6000 sec for both the blue and red arm. 
We used IRAF for a standard reduction of the data (bias 
subtraction, flatfielding, wavelength calibration, cosmic-ray 
removal, background subtraction and tilt removal). Arc ex- 
posures for wavelength calibration were taken before and 
after the observations of B2 0722+30 at approximately the 
same position. We used a 60 sec observation of the telluric 
standard star HD58336 to correct for atmospheric absorp- 
tion bands in the red part of the spectrum as well as possible. 
The resulting spectra are aligned within one pixel in the spa- 
tial direction. For the flux calibration we used five standard 
calibration stars (D191-B2B, Feige 67, PG0216+032, Feige 



24 and HD19445). The accuracy of the relative flux calibra- 
tion is within 6% (only in the very red part of the spectrum, 
beyond 7400 A, the flux calibration errors are up to 10%). 
Due to slit-losses when observing the flux-calibration stars 
we could not normalise the spectrum better than within a 
factor of 2. We note that this does not affect the stellar 
population analysis of Sect. 13.21 (because the relative flux 
calibration is unaffected) , but it does introduce a factor of 2 
uncertainty in the total stellar mass estimates derived from 
the spectra fSect. [3~2.3|l . Subsequently, we used the Starlink 
package FIGARO to correct the spectra for galactic extinc- 
tion (assuming E(B - V) = 0.0590 and to de-redshift the 
spectra to rest-wavelengths. For the analysis of the spectra 
we used the Starlink package DIPSO for emission-line fitting 
and the Interactive Data Language (IDL) for the stellar pop- 
ulation modelling. 

2. 3. 1 Stellar population synthesis method 

To investigate in detail the stellar populations in 
B2 0722+30, we model the continuum Spectral Energy Dis- 
tribution (SED) of the optical spectra in various regions, 
taking into account both stellar and AGN-related compo- 
nents. Subsequently, we will make a more detailed compar- 
ison between the data and the models by investigating the 
age-sensitive Call K and Balmer absorption lines in order 
to constrai n our result better. This techn i que has been used 
befor e (see iTadhunter et"afl l200l 120051 : IWills et al ] I20021 . 
l2004h and is preferred to using absorption line indices at 
face-value, because most of the age sensitive diagnostic ab- 
sorption lines could be affected by emission-line contamina- 
tion due to the activity in B2 0722+30. 

For mo delling the observed spectr a, stellar population 
models from lBruzual fc Chariot! |2003) are used. These are 
based on Salpeter initial mass function (IMF) and solar 
metallicity, instantaneous starbursts. We use a x 2 minimi- 
sation technique to fit combinations of an unreddened 12.5 
Gyr old stellar population (OSP), a young stellar population 
(YSP) and a power-law component of the form F\ tx \ a 
to the observed spectra (see Sect. 13.21 for more details). 
YSP template spectra across different ages and with various 
amounts of reddening are used for this. The latter is nec- 
essary, since we were not able to determine the reddening 
a priori using Hq/H/3 in B2 0722+30, because H/3 absorp- 
tion appears to be significant in this system, which dilutes 
the H/3 emission line (and apart fr om Ha, t he ot her Balmer 
lines are too weak to be used). The Scaton (1979) reddening 
law is used to apply the reddening to the template spectra. 
The final library of YSP template spectra comprises YSP 
models with age 0.01 - 5.0 Gyr and with E(B - V) = 0.0 

- 1.6. We compare the total flux of the combined OSP plus 
YSP (plus possible power-law component) model with the 
observed flux in wavelength-bins across the spectrum. In or- 
der to scale the template models to the observed spectra, 
a normalising bin was chosen in the wavelength range 4720 

- 4820 A. For the x 2 -ntting we assume an error of ±6% 
in each wavelength bin, in agreement with the uncertainty 
in the flux calibration. Note that, since the flux calibration 

2 Based on results from the NASA/IPAC Extragalactic Database 
(NED) 
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errors are not likely to be independent between the various 
wavelength bins, we can merely use the reduced x 2 values as 
an indication of the region of parameter space for which the 
modelling provides good results, rather than derive accurate 
statistical properties of the fitting procedure itself. For that, 
we also need to inspect the model-fit to our spectra visually. 
The Ca II K (A 3933 A) and higher order Balmer absorption 
lines in the blue part of the spectrum, as well as the G-band 
feature (A 4305 A), are excellent to use for this. 

As investigated bv lDickson et al" I <| 19951 . and references 
therein), spectra of active galaxies may be diluted by UV- 
excess due to nebular continuum and AGN related effects 
in the region sh ortwa rd of 4000 A ( see also discussion in 
Tadhunter et alj|2002f ). In a previous paper (jEmonts et al.l 
20061 ) we investigated in detail the effects of nebular contin- 
uum in the spectral light of radio galaxy B2 0648+27, which 
was observed during the same observing run and under the 
same conditions as B2 0722+30. We concluded that for the 
nuclear spectrum of B2 0648+27 the most extreme cases - 
maximum nebular continuum subtraction and no nebular 
continuum subtraction - did not significantly alter our main 
results of the SED modelling. The emission-line equivalent 
width provi des an indication for the contribution of nebular 
continuum (|Dickson et al. 19951) as well as s cattered AGN 
components (see e.g. iTadhunter et al.ll2002l ). B2 0722+30 
does not contain emission lines that are significantly stronger 
than the emission lines in B2 0648+27. We therefore do not 
expect that nebular continuum will have a significant con- 
tribution to the light in the spectra of B2 0722+30, hence 
no correction for the nebular continuum is applied. 



disc is centred on 5675 km s . This is in agreement with 
the central part of the absorption as well as with redshift 
det erminations f r om C O data by iMazzarella et all (Il993l ) 
and lEvans etall (|2005l ). 

Figure H shows a plot of the H a +[N II] AA6548,6583 
emission lines that we detect along the major axis of 
B2 0722+30. Also here, regular rotation is clearly visible. 
We trace the emission lines out to about 5.4 kpc on either 
side of the nucleus. The optical emission lines cover total 
velocity range of 322 ±11 km s _1 and are centred on a 
redshift velocity of 5666 ± 60 km s -1 (which corresponds to 
z — 0.0189 ± 0.0002). We assume this is the systemic ve- 
locity of B2 0722+30, which is in good agreement with the 
central velocity of the H I and CO gas described above. 

At the northwestern tip of the disc there appears to be 
a tail of H I gas stretching to the north, in the direction of an 
optical counterpart (marked as "XI" in Figs. [1] and The 
radial profile of this optical counterpart (deduced with the 
interactive graphics anaysis tools in IRAF) is significantly 
broader than that of the stars in the field. We therefore argue 
that this counterpart likely represents a small galaxy rather 
than a foreground star. If it is a galaxy, it can be either a true 
companion of B2 0722+30 or a galaxy in projection. Up to 
48 kpc further out to the north (around dec=29°59') there 
appears to be more faint HI emission, although this has to 
be confirmed with additional data. Another, even fainter, 
peculiar optical feature is present about 6 kpc southwest of 
B2 0722+30. This feature is only visible in the deep optical 
image of Fig.[T]("X2"). The nature of this very faint optical 
feature is not clear. 



3 RESULTS 

Figure[T]shows the deep B-band image of B2 0722+30 and its 
environment. B2 0722+30 (galaxy 'A') is clea rly dominated 
by a s tellar disc. The optical HST imaging bv lCapetti et al.l 
( 2000) shows in more detail the bulge component and dust 
features. Several other interesting galaxies and optical fea- 
tures are visible in the environment of B2 0722+30 in Fig. [T] 
They have a striking connection with the neutral hydrogen 
gas in the region around B2 0722+30 and will therefore be 
discussed in detail in the next Section. 



3.1 Neutral and ionised gas 

3.1.1 Host galaxy 

Figure [2] shows a total intensity map of the HI gas detected 
in B2 0722+30 and its environment, while Fig. [3] shows the 
channel maps of the HI data. The HI gas in B2 0722+30 
is distributed in a rotating disc that follows the optical disc 
of the host galaxy out to a radius of about 8 kpc from the 
centre. The HI mass detected in emission in this disc is 
Mm = 2.3 x 1O 8 M0 and the average surface density is about 
4.1 Mq pc -2 . The total H I mass in the disc will be somewhat 
larger than this, because the part of the HI gas that lies in 
front of the radio continuum source is detected in absorption 
and not included in this HI mass estimate. The HI disc 
seems to be in regular rotation and covers a velocity range 
of 450 km s~\ which is equal to the FWHM of the CO 
emission-line profile detected bv lEvans et al. (2005). The HI 



3.1.2 Environment 

The environment of B2 0722+30 shows a number of lumi- 
nous galaxies that appear to be in interaction. The most in- 
teresting is the galaxy pair UGC 3841 east of B2 0722+30. 
Our deep optical B-band image (Fig. [TJ suggests that the 
western galaxy of the UGC 3841 pair (at a distance of 36 
kpc from B2 0722+30) is most likely a dust-lane ellipti- 
cal. This galaxy contains a weak, unresolved radio contin- 
uum source, MRK 1201, which h as Seyfert 2 characteristics 
ijDessauges-Zavadskv et al]|2000h . From our 1.4 GHz contin- 
uum data we derive a flux of Fxa GHz = 3.1 mjy beam - 
(Pi. 4 GHz = 2.3 x 10 21 W Hz -1 ) for this source. For clar- 
ification, we refer to this western galaxy of the UGC 3841 
pair as MRK 1201. The eastern galaxy of this pair (which 
we will refer to as UGC 3841a) is an edge-on, distorted 
and elongated disc galaxy. The western part of the disc of 
UGC 3841a extends past MRK 1201 and points in the direc- 
tion of B2 0722+30. MRK 1201 and UGC 3841a appear to 
be in heavy interaction, as can be seen from the kinematics 
of the H I gas in Fig. [3] The H I gas follows the distorted disc 
of UGC 3841a, except in the velocity range 5629 ^ v ^ 5800 
km s _1 , where it is spread around MRK 1201. An HI bridge, 
about 22 kpc in length, seems to stretch from the HI dis- 
tribution around MRK 1201 in the direction of B2 0722+30 
(5650 ^ v ^ 5714 km s _1 in Fig. [3}. The total HI mass 
detected in the MRK 1201/UGC 3841a pair is 3.4 x 10 9 M Q . 

To the north, at a distance of 33 kpc from UGC 3841a, is 
another Hi-rich galaxy pair. To our knowledge, this galaxy 
pair has not been classified previously. We refer to these 
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Figure 1. Deep optical B-band image of B2 0722+30 and its environment. The high-contrast plot (left) marks various galaxies and 
features of interest that are plotted in low contrast on the right and which are discussed in more detail in the text. 
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Figure 2. Left; total intensity map of HI emission in B2 0722+30 and its environment (contours) overlaid on to the optical MDM image 
(grey-scale). Contour levels: 0.63, 1.2, 1.8, 2.7, 3.3, 4.3, 5.3, 6.3, 7.7 xlO 20 cm -2 . The zoom-in shows B2 0722+30 in detail (for clarity, 
the part of the gas-disc that is seen in absorption against the radio source is not shown). Right: Position- velocity plot along the major 
axis of the HI disc, as indicated in the zoom-in on the left. Contour levels: -0.5 -1.4 -2.4 -3.4 -4.4 (grey); 0.5 0.7 0.9 1.1 1.3 (black) mjy 
beam -1 . The profile of the HI absorption against the radio continuum is also show. The bar indicates the systemic velocity. 



galaxies as J072549+295830 (east) and J072548+295828 
(west). This galaxy pair contains 4.2 x 10 8 M Q of HI. 



3.2 Stellar populations 

In order to study the stellar populations throughout 
B2 0722+30, we extracted spectra at various apertures along 
the optical emission-line disc. Figure [4] shows the location 
and Table [1] gives the details of these apertures. In Table 
[1] columns 2 and 3 give the size of the aperture (in arcsec 
and kpc), while columns 4 and 5 give the distance of the 
centre of the apertures from the nucleus (assumed to be the 
brightest peak in the stellar continuum of the spectrum). 



Table 1. Apertures 





Width 


Width 




^nuc 




(arcsec) 


(kpc) 


(arcsec) 


(kpc) 


NUC 


2.0 


0.8 






SE1 


3.2 


1.2 


1.2 


1.6 


SE2 


6.0 


2.3 


9.2 


3.6 


NW 


6.0 


2.3 


7.2 


2.8 



The method we used for the stellar population analysis in 
the various regions is described in Sect. l2.3.TI 
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Figure 3. Channel maps of the environment of B2 0722+30 from the smoothed data cube. The HI bridge between MRK 1201 and 
B2 0722+30 is visible in the channels with 5714 v ^ 5650 km s" 1 . Contour levels: -0.7, -1.7, -3.0, -4.5 (grey); 0.7, 0.9, 1.1, 1.4, 1.8, 2.3, 
2.9, 3.7 (black) mjy beam -1 . The first panel shows the total intensity image of Fig.[2](i.e. higher resolution than the channel maps). 



3.2.1 Nucleus 

First we model the spectrum of the nuclear aperture of 
B2 0722+30. A single unreddened 12.5 Gyr OSP does not 
provide an accurate fit to the overall SED in the nuclear 
spectrum, i.e. the residuals of the model-fit to the spectrum 
are much larger than the flux calibration errors (Sect. [2~3)l . 
When including a power-law component in addition to the 
12.5 Gyr OSP, a good fit can be obtained only for a power- 
law component with a — 2.8, which contributes to about 
50% of the light in the red part of the spectrum. A power- 
law with such a slope and high luminosity is not likely the 
result of an optical AGN in B2 0722+30, because there is no 
cl ear evidence for a br ight point-source in the HST image 
of ICapetti et al.l (|2000h and because (as we will see in Sect. 
I3.3|l the nuclear emission lines are characteristic of a Low 
Ionization Nuclear Emission-line Region (LINER) rather 
than an AGN power-law model. To investigate whether the 
power-law component could represent an additional redden- 
ing component that we did not take into account, we tried 



to fit a 12.5 Gyr OSP that we reddened by various amounts, 
but this did not produce a good fit. 

Since no good or physically acceptable fit of a single 
12.5 Gyr OSP or OSP + power-law component could be 
obtained for the nuclear region of B2 0722+30, we subse- 
quently modelled the nuclear spectrum with a combination 
of a 12.5 Gyr OSP and a (reddened) YSP. Figure [5] shows 
the results of adding a YSP. The solution of fitting an OSP 
+ YSP clearly converges to two different models: including a 
YSP of age 0.05 Gyr, E(B-V) = 1.2 and light percentage = 
30% results in reduced % 2 = 0.27, the least x 2 fit across our 
range of parameter space; including a YSP of age 2.0 Gyr, 
E(B — V)— 0.4 and light percentage = 89% gives a least \ 2 
solution of 0.52 in this part of the parameter space. These 
two model fits are shown in Fig. [5] In an attempt to break 
this degeneracy, we look in detail at the age sensitive absorp- 
tion features around 4000 A. Figure [7] (top) shows that both 
models provide a reasonable fit to the age sensitive diagnos- 
tic Balmer and Call K(A3933)/CaII H(A3968) absorption 
lines and the G-band feature around A 4305 A. Overall, the 
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Figure 4. 2D-spectrum of B2 0722+30 around 
Ha+[N II] AA6548,6583 (contour levels are between 10 and 
90% of the maximum in steps of 5%). On the right the various 
apertures that we used for our SED modelling are indicated. 



0.05 Gyr YSP seems to provide the best fit to most of the 
features, except for Call K. Narrow Call K absorption by 
the interstellar medium (ISM) might affect t he line-strength 
in th e core of the stellar Call K line (e.g. iPasquini et al.l 
2000), resulting in a deeper Call K absorption feature in 
the observed spectra compared with our model. The fit of 
an intermediate age stellar population with moderate red- 
dening (age 2.0 Gyr and E(B — V) = 0.4) provides a better 
relative fit to both calcium lines (Figure [7]- bottom), but 
appears to slightly over predict the depth of the lines. We 
note, however, that part of this could be due to a slight 
difference in resolution between the real and the template 
spectra. We also tried to fit a two-component combination 
of a 2 Gyr OSP plus 0.05 Gyr YSP to the nuclear region, 
but due to the dominating contribution of the 2.0 Gyr stel- 
lar population in that case (>85% of the total stellar light) 
this does not provide accurate information about a possible 
YSP component. Including an additional power- law compo- 
nent to our model fits did not significantly improve or alter 
our results. 

We therefore argue that our used method cannot clearly 
distinguish between a dominating 2.0 Gyr stellar population 
or a combination of a 12.5 Gyr OSP and a heavily reddened 
0.05 Gyr YSP in the nuclear region of B2 0722+30 (which 
will be dominated by the bulge component). Further limita- 
tions of our used method are discussed in Sect. 13.231 



3.2.2 Off-nuclear regions 

The spectra in the off- nuclear regions (SE1, SE2 and NW) 
cannot be fitted accurately with a single unreddened 12.5 
Gyr OSP and an AGN power-law component is not likely to 



be important far from the nucleus (which did not need the 
inclusion of a power-law component anyway). 

Figure [5] shows the results of adding a YSP in fitting 
the spectrum of region SE1. As can be seen from 'case A', 
when adding a YSP of age ~ 0.05 -0.1 Gyr and E(B - V) « 
0.4 — 0.7 to the 12.5 Gyr OSP, we get a reasonable fit to 
the SED in the region SE1. However, when looking in detail 
at the model fit of case A (12.5 Gyr OSP + YSP), we see 
that the fit is not very accurate to the age-sensitive Ca II K 
line. Although narrow Ca II K absorption by the ISM might 
affect th e line strength in th e core of the stellar Call K 
line (e.g. IPasquini et al.ll2000l ). the asymmetry between the 
observed spectrum and our model fit for the Call K line 
in the off-nuclear region SE1 is much larger than for the 
nuclear region. We do not expect that this can be explained 
solely by Call K absorption by the ISM. This questions the 
validity of our SED fitting results in region SE1. However, 
we saw that the nuclear region may be dominated by an 
intermediate age stellar population of 2.0 Gyr instead of a 
12.5 Gyr OSP. We therefore also modelled our off-nuclear 
spectrum in region SE1 with an unreddened 2.0 Gyr "OSP" 
combined with a YSP ('case B' in Fig. [8}. Substituting an 
underlying 12.5 Gyr OSP for a 2.0 Gyr OSP gives a much 
better fit to the Call K line. At the same time, the least 
\ 2 solution for the overall SED fit does not change much. 
For the case B modelling the YSP has an age of about 0.05 
(NW) - 0.1 (SE1 and SE2) Gyr and contributes to about 
35-55% of the total stellar light. The reddening is slightly 
higher than for the case A modelling. 

The story for regions SE2 and NW is very similar to 
that of region SE1. Also for these regions, the combination 
of a 12.5 Gyr OSP + YSP gives an acceptable least \ 2 so ~ 
lution, but again severely under-predicts the strength of the 
Call K line. Substituting the 12.5 Gyr OSP for a 2.0 Gyr 
OSP provides a good fit to the age-sensitive Call K-line, 
as well as the Balmer absorption lines and G-band feature, 
without changing the age and light contribution of the YSP 
dramatically. The best x 2 solution and fit to the spectra in 
regions SE2 and NW are shown in Fig. 

Therefore - although there is a substantial uncertainty 
in the properties of the OSP across B2 0722+30 as well as 
the stellar content of the bulge region - the presence of a 
YSP across the disc of B2 0722+30 is fairly well established 
with our stellar population method. Table [5] summarises 
the best-fit results from our stellar population analysis of 
B2 0722+30. The light percentage coming from the YSP is 
roughly 40-50% for all the three off-nuclear regions. This 
strongly suggests that the stellar population is fairly uni- 
form across the disc of B2 0722+30. It is interesting to note 
that the reddening of the YSP, E(B — V), is highest in the 
inner region SE1 {E(B — V) ~ 1.0) and gets lower for the 
outer regions (E(B - V) » 0.7 in NW and « 0.5 in SE2). 

3.2.3 Stellar masses 

The masses of the OSP and YSP in the off-nuclear apertures 
are also given in Table[5] The masses are calculated from the 
flux of the YSP and OSP at 4770 A (the central A of our 
normalising bin). The YSP and OSP fluxes (taking into ac- 
count reddening; see Table [ ^jl are scaled to stellar mass es us- 
ing the template spectra of Bruzual fc Chariot! ([2003) . The 
uncertainty in the stellar mass estimates is based on uncer- 
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Figure 6. B2 0722+30 NUC - Least \ model fits to the nuclear spectrum of B2 0722+30. The dotted line is the observed spectrum, 
the black solid line is the fitted model and the grey solid lines are the separate components of the fitted model. The histogram at the 
bottom of the plots shows the residuals of the model fit to the observed spectrum. Top: combination of a 12.5 Gyr old stellar population 
(OSP) and a young stellar population (YSP) of age = 0.05 Gyr and E(B — V) = 1.2. Bottom: combination of a 12.5 Gyr old stellar 
population (OSP) and an intermediate age stellar population of age = 2.0 Gyr and E(B — V) = 0.4. 



tainties from the SED modelling and does not include the 
uncertainty of a factor of 2 due to slit-losses during the ob- 
servations (see Sect. 12.31 for details). 

3.2-4 Infrared properties 

B2 0722+30 has a far-infrared (FIR) luminosity of Lpir = 
4.7 x 1O 1O L l|Mazzarella et al.lll993l . corrected for H = 71 



km s _1 Mpc~ ), similar to that of star forming Sb and Sc 
galaxies l|Rieke fc Lebofskvlll986l ). This is a relatively high 
FIR luminosity, given that - as we will see in Sect. 14.11 - 
B2 0722+30 is relatively small compared with for exam- 
ple the Milky Way galaxy. In order to investigate whether 
dust-heating by the YSP in B2 0722+30 could be respon- 
sible for this high FIR luminosity, we calculate the total 
bolometric luminosity (Lboi) of the 0.05-0.1 Gyr YSP across 
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Figure 7. B2 0722+30 NUC - Least \ 2 model fits to the nuclear 
spectrum of B2 0722+30 in the region around 4000 A. The dotted 
line is the observed spectrum, the black solid line the fitted model. 
Top: combination of a 12.5 Gyr old stellar population (OSP) and 
a young stellar population (YSP) of age = 0.05 Gyr and E(B — V) 
= 1.2. Bottom: combination of a 12.5 Gyr old stellar population 
(OSP) and an intermediate age stellar population of age = 2.0 
Gyr and E(B - V) = 0.4. 



B2 0722+30. Lboi is estimated from the integrated lumi- 
nosity of the unreddened YSP measured from our template 
spectra (across the wavelength range - 30,000 A), scaled 
to the total mass of the YSP given in Table [2] Given that 
the reddening is E(B — V) = 0.5 — 1.0 across the disc, it is 
likely that most (> 85%) of the total bolometric luminosity 
of the young stars is absorbed by dust. Assuming that most 
of this dust-absorbed light is re-radiated into the FIR, Lboi 
should therefore provide a reasonable estimate for the FIR 
luminosity related to the YSP in B2 0722+30. 

We estimate that L bo i of the YSP in B2 0722+30 is be- 
tween 0.3 x 10 10 L Q and 1.5 x 10 10 L Q , depending on the age 
of the YSP varying between 0.05 and 0.1 Gyr and whether 
or not a 0.05 Gyr YSP in the nuclear region is taken into 
account (see Table [2}. We argue that these estimates are 
likely lower limits, given that our apertures do not cover the 
entire optical extent of the host galaxy. In addition, stars 
which formed less than 0.05 — 0.1 Gyr ago (which are not 
identified with our method for modelling the stellar popula- 
tion) provide the highest relative contribution to the YSP's 
bolometric luminosity. We note, however, that this estimate 
does not include the factor 2 uncertainty in absolute flux 
calibration of our spectra (see Sect. 12. 3p 

In order to circumvent the limited slit coverage and un- 
certainty in absolute flux calibration and to obtain a bet- 
ter estimate of the total bolometric luminosity of the en- 
tire YSP in B2 0722+30, we apply the average values of 
our SE P modelling to the total V-band lum inosity of the 
system. iGonzalez-Serrano fc Carballol l|2000l ) measured the 
optical V-band magnitude of the host galaxy B2 0722+30 
to be mv = 15.04, which corresponds to Lv = 5.0 x lO 9 Z/0. 
Assuming that 43% of the V-band light across B2 0722+30 



comes from a 0.05-0.1 Gyr YSP with reddening E(B — V) = 
0.73, we estimate that Lboi = 3.5 — 5.0 x 1O 1O L0, which is 
in excellent agreement with Lpir measured from the IRAS 
data. We therefore argue that the total Lboi of the YSP 
appears to be roughly consistent with Lpir in B2 0722+30 
and that dust heating by the YSP in B2 0722+30 is likely 
responsible for at least a substantial fraction of the relatively 
high FIR luminosity of the host galaxy. 

3.2.5 Uncertainties in our method 

Of course there is a degree of uncertainty introduced by ob- 
servational errors (such as flux calibration), the uniqueness 
of the model solutions (see Sect. l3.2"Tl and [3,2.2[) and by the 
assumed parameters of the synthesis models that we used in 
our fitting procedure (duration of the starburst, the assumed 
Salpeter IMF and metallicity). 

Uncertainties in the flux calibration could be the rea- 
son that in particular in the red part of the spectrum the 
fit is not perfect. However, we do not expect that this will 
significantly change our results, since our visual inspection 
of the detailed fits in the blue part of the spectrum (around 
4000 A) is in excellent agreement with the reduced % 2 re- 
sults. There is also the selection effect that the apertures 
that we chose for our study cover regions of a kpc or more 
in various parts of the galaxy's disc. Locally, the stellar pop- 
ulations and dust properties could vary from place to place 
within these apertures. Nevertheless, the similarity of the 
results for the various off-nuclear regions reassures us that 
outside the nucleus the average values for the stellar ages 
and masses are fairly uniform. All four apertures combined 
do not completely cover the total extent of the galaxy, there- 
fore the total mass of the stellar populations in B2 0722+30 
is somewhat larger than the sum of the stellar masses across 
the separate apertures. 

For a detailed discussion of the uncertainties in the 
assumed parameters o f the synthesis models we refer to 
iTadhunter et all (|2002h . As discussed there, the major un- 
certainty in the model parameters is the assumed shape of 
the IMF. This may lead to uncertainties of a factor of 2- 
3 in the total mass estimates in Table [2] (on top of the 
factor of 2 due to observational inaccuracies; Sect. 12. 3p . 
Regarding the metallicity, we only used template spectra 
with solar-abundances t o fit ou r data. Following for example 
Sarichcz-Blazq uez et al.l (|2006l ). uncertainty in the metallic- 
ity could cause a serious scatter in the derived Single Stel- 
lar Population (SSP)-equivalent age of a stellar population. 
However, we do not expect that metallicity effects signifi- 
cantly influence the age and mass estimates that we derive 
for the YSP detected in th e disc of B2 722+ 30, where one 
expects solar metallicities. ISarzi et"aL (2005) investigated 
the effect of metallicity on an SED fitting technique for the 
bulge-regions of nearby galaxies. They show that substitut- 
ing solar for 2.5 x solar metallicity YSP template spectra 
does not change the results for age and light percentage of 
the YSP dramatically. The most significant effect on the 
YSP is that super-solar metallicities can mimic additional 
reddening, in which case our solar-abundance models might 
overestimate the derived E(B — V). Metallicity effects could, 
however, seriously affect the results for the older 2.0 Gyr 
stellar population that we fit throughout B2 0722+30 and 
which may dominate the stellar mass in the bulge region. It 
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Figure 8. B2 0722+30 SE1 - Model 'A' includes an unreddened OSP of 12.5 Gyr; model 'B' an unreddened OSP of 2.0 Gyr. Top left: 
least x 2 results of the SED modelling for different age and E(B — V) of the YSP. Contour levels are the same as in Fig. \E\ Top right: 
Detailed least \ 2 model fit in the region around 4000 A. The dotted line is the observed spectrum, the black solid line the fitted model. 
Bottom: Least x 2 model fit to the entire SED. The dotted line is the observed spectrum, the black solid line the fitted model and the 
grey solid lines the separate components of the fitted model. The histogram at the bottom of the plots shows the residuals of the model 
fit to the observed spectrum. 



Table 2. Best results of the least \' 2 model fitting 



Aperture 


Age OSP 


Age YSP 


E(B-V) 


Light % 


Mass % 


Mass (M ) 


Mass (Mq 




(Gyr) 


(Gyr) 


YSP 


YSP 


YSP 


YSP 


OSP 


NUC (1) 


12.5 


0.05 


1.2 


30 


22 


4.1 x 10 8 


3.4 x 10 9 


NUC (2) 


12.5 


2.0 


0.1 


89 


83 


2.4 x 10 9 


5.3 x 10 8 


SE1 


2.0 


0.1 


1.0 


37 


61 


2.1 x 10 8 


2.7 x 10 8 


SE2 


2.0 


0.1 


0.5 


48 


29 


8.4 x 10 7 


2.8 x 10 s 


NW 


2.0 


0.05 


0.7 


11 


.32 


8.9 x 10 7 


3.0 x 10 s 



Notes — The light% is the relative flux contribution of the YSP in the normalising bin (4720 - 4820 A). The mass of the stellar 
populations in the various apertures (last two columns) is derived from the stellar population's light% of the observed flux at 4800 A - 
corrected for reddening - and scaled to the unreddened model template spectra. 



is known that ther e is a broad range in metallicities across 
galaxy bulges (e.g. IWvse et al.lll997T ) and that metallicities 
in the central bulge can be significantly higher than in the 
disc. Moreover, for old and intermediate age stellar popu- 
lations such as the older stellar population in B2 0722+30 
an age-metallicity degeneracy is likely to be much more sig- 
nificant than for young stellar populations. Therefore, our 
results for the 2.0 Gyr stellar population in B2 0722+30 
should be taken with care and we warn the reader that the 
estimated age of this older stellar population should not be 
take at face-value. Finally, our used models are based on 



fitting a single, instantaneous starburst event on top of an 
old stellar population. More subtle, previous starbursts or 
longer periods of continuous star formation cannot be de- 
tected with our method. It is very likely that the YSP de- 
tected in the disc of B2 0722+30 is evidence for ongoing star 
formation over a substantial period, as is normal in spiral 
galaxies. 
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Figure 9. B2 0722+30 SE2 and NW Best fit of a 2.0 Gyr OSP + reddened YSP to the spectra in regions SE2 and NW. Top: Least \ 2 
model fit to the entire SED. The dotted line is the observed spectrum, the black solid line the fitted model and the grey solid lines the 
separate components of the fitted model. The histogram at the bottom of the plots shows the residuals of the model fit to the observed 
spectrum. Bottom: Detailed least y 2 model fit in the region around 4000 A. The dotted line is the observed spectrum, the black solid 
line the fitted model. 



3.2.6 Stellar populations in B2 0722+30 (summary) 

The analysis of the spectra of B2 0722+30 shows that a 0.05 
- 0.1 Gyr young stellar population contributes to about half 
of the stellar light in the disc of B2 0722+30. This YSP most 
likely traces ongoing star formation (as we will see in Sect. 
13.31 in the disc of B2 0722+30 the emission-line ratios re- 
semble those of star- forming regions). We argue that there is 
a significant uncertainty only in the total mass and perhaps 
in the reddening of this young stellar population. The stellar 
content of the bulge region as well as the properties of the 
underlying old stellar population throughout the galaxy are 
much more uncertain and could not be properly constrained 
with our used SED modelling technique. 



3.3 Emission-line diagnostics 

In order to investigate the optical properties of the AGN 
and to get an idea about the ionisation-mechanism of the 
emission-line gas in B2 0722+30, we measured the emission- 
line ratios in the nuclear as well as the off- nuclear regions. 
Emission-line infilling by stellar-related absorption features 
(in particular in the hydrogen lines) could be significant, 
therefore we subtracted the best-fit model from our observed 
spectra to eliminate the underlying stellar continuum and 
absorption as well as possible. Table [3] gives the emission- 
line ratios in the various apertures. The last column gives 
the corresponding charact eristic spectral classification b ased 
on diagnostic d iagrams of IVeilleux fc Qsterbrockl (I1987T ) and 
iBaldwin et al.l (1 198 lT) . In p articu lar the diagnostic diagrams 
of IVeilleux fc Qsterbrockl l|l987l ) are based on the ratios of 



extinction of E(B - V) = 0.059; see Sect. [23]) is L [0 iii] = 
1.3 x 10 39 erg s" 1 . 

The nuclear region shows the emission-line character- 
istics of a Low Ionisation Nuclear Emission-line Region 
(LINER). The LINER classification for B2 0722+30 can 
clearly be distinguished from a Seyfert-type nucleus i n the 
diagnostic dia g rams of IVeilleux fc Qsterbrockl l|l987l ) and 
IBaldwin et al.l (| 198 il l. We note, however, that our nuclear 
aperture covers a region of about 0.8 kpc, which is sig- 
nificantly larger than the central black-hole region. Never- 
theless, a bright Seyfert nucleus would likely dominate the 
emission-line flux in this nuclear aperture and no indications 
of such a strong optical AGN are present in our spectra. 
The emission-line spectrum of the nuclear region therefore 
indicates that the optical AGN has the characteristics of a 
low-excitation AGN, as is common for FR-I radio galaxies. 
We will discuss this further in Sect. 14.21 

The off nuclear regions in the disc of B2 0722+30 show 
emission-line characteristics that resemble those of HII re- 
gions. This strengthens the idea that the YSP detected 
across the disc of B2 0722+30 represents ongoing star for- 
mation. 



emission lines that are close together in order to minimize 
the effects of reddening. The [OIII]A5007 luminosity that 
we measure in the nuclear aperture (corrected for a galactic 



4 DISCUSSION 

In this Section, we discuss in detail the properties of the host 
galaxy and AGN in B2 0722+30. We compare the charac- 
teristics of B2 0722+30 with two other disc-dominated radio 
galaxies in order to investigate in more detail the circum- 
stances under which a double-lobed radio source can occur 
in these disc-galaxies. 
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Table 3. Line ratios 



Aperture 


[O III] 

H f> 


[OH] 
H P 


[OH] 
[OIIIJ 


[Nil] 
H a 


[SII] 
H a 


[O I] 
H a 


type 


NUC 


1.5 


3.3 


2.2 


1.2 


0.8 


0.14 


LINER 


SE1 


0.3 


0.9 


2.(5 


0.5 


0.3 


0.04 


HII region 


SE2 


0.4 


1.3 


3.2 


0.5 


0.3 


0.06 


HII / LINER 


NW 


0.8 


1.7 


2.1 


0.6 


0.4 


0.05 


H II / LINER 



Notes - Line-ratios are given for the emission lines H/3A4861, HaA6563, [OIII]A5007, [OII]A3727, [NII]A6583, [S II] A6716+A6731 and 
[OI1A6300. The last row gives the characteristic type of emission-line object that corresponds to the line-ratios (as given by 
iBaldwin et~allll98ll; IVeilleux fc Osterbrocklll987h . 



4.1 Disc-dominated host galaxy 

The optical morphology, gas content/distribution and stel- 
lar population analysis described in this paper clearly show 
that B2 0722+30 is a galaxy with a prominent star-forming 
disc. The optical morpholo gy of B2 0722+30 is m ost clearly 
seen in the HST image of iCapetti ei~al] (|2000h . who clas- 
sify the galaxy as a spiral. However, given the large inclina- 
tion of the galaxy, direct evidence for a spiral morphology 
in the form of prominent spiral arms is lacking and hence 
the classification of B2 0722+30 is often ambiguous in the 
literature (spiral or SO). From V-ba nd surface photometry, 
iGonzalez-Serrano fc Carballol (I2000T I classify B2 0722+30 as 
a spiral or SO galaxy and conclude that its radial-light pro- 
file has discy isophotes that imply the presence of either a 
disc, arms or a bar. 

Figure [TU] shows in more detail our deep optical B-band 
image of B2 0722+30. There is clear morphological structure 
visible a long the stellar disc . A comparison with the HST 
image of iCapetti et alj f2000) shows that at least part of this 
structure is caused by the presence of large amounts of dust 
across the disc, although the exact extent of the dust remains 
somewhat ambiguous from our data (as well as from the 
limited field-of-view of the HST image). The presence of a 
gas-rich disc with large-scale dust and active star formation 
fits the classical definition of a late-type system, although we 
note that similar properties have recently been observed for 
NGC 612 (|Holt et alj|2007l ; lEmonts et al.ll2008l ). which was 
classified as an archetypal S O from its optical morphology 
l|Veron-Cettv fc Veror]|200ll ). Figure \W\ also shows that the 
bulge region of B2 0722+30 (marked by visually comparing 
our B-band image with the HST image of Capetti et al.) is 
modest compared to the total extent of the disc. We derive 
an estimated bulge-to-disc flux ratio of B/D ~ 0.2 from our 
B-band image (which is likely an upper limit, given the large 
inclination and dust-obscuration of the disc). This B/D is in 
agreement with what is generally found for spiral galaxies 
(although the average B/D for SOs is only slightly high er 
than our estimated upper limit; iGraham fc Worlevl l2008). 



From the observational evidence, we therefore conclude 
that B2 0722+30 is a gas-rich disc galaxy. Although its char- 
acteristics appear to fit the description of a spiral galaxy, the 
exact morphological classification is not expected to influ- 
ence the main conclusions of this paper, given that for for 
disc galaxies in general (both spirals and SOs) it is extremely 
rare to host a classical radio source (see Sect. [1] for details). 

Our HI and optical emission-line data show that 
B2 0722+30 contains a regularly rotating disc of both cool 
neutral hydrogen gas and ionised emission-line gas. The 




Figure 10. Deep optical B-band image of B2 0722+30. Contours 
are in arbitrary count-units and range from 5 to 900 in steps of 
100. The innermost contour marks the bulge region, based on 
visual comparison with the HST image of Capetti ct al. (2000). 



total mass of neutral hydrogen gas in the disc (A/hi = 
2.3 x 10 s Mq, or a somewhat higher value when taking into 
account the significant amounts of HI gas detected in ab- 
sorption against the radio continuum) is a factor of five 
lower than the mass of molecular hydrogen in B2 0722+30 
(M H 2 = 1. 2 x 10 9 M ) 7), as deri ved from single-dish CO obser- 
vations bv lEvans et al.l [20051 ). and about an order of mag- 
nitude lower than the H I mass of our Milky W ay galaxy 
(Mm mw = 5 x 10 9 M Q ; iHenderson et alj|l982h . Also, the 
total extent of the HI disc in B2 0722+30 (~ 15 kpc) is less 
than half the diameter of the Milky Way HI disc. As men- 
tioned in Sect. 13.2.41 t he V-band luminosity of B2 0722+3 
is L v = 5.0 x 10 9 L Q l|Gonzalez-Serrano fc Carballol l200Ch ■ 
Therefore, Mm/Ly = 0.05. Albeit at the low end, this 
is consistent for w hat is observed in spiral galaxies (e.g 
iKnapp et all 1 19891 . please note that our derived value of 
Mm/Ly is somewhat larger when taking into account the 
H I gas detected in absorption) . The average surface density 
of the disc is about 4.1 Mq pc~ 2 , which is close to the crit- 
ical gas surface density fo r star formation in galaxy discs 
ijMartin fc Kennicuttll200ll : Ivan der Hulst et al.lll993l ). The 
YSP across the disc and the emission-line ratios character- 
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istic of H II regions in the ionised gas disc indicate that star 
formation is happening across the disc of B2 0722+30. 

Given the regular rotation of the HI gas in the disc of 
B2 0722+30, we can make an estimate for the total mass 
enclosed by the system by assuming that the underlying 
dark matter halo has a spherical distribution. In that case, 



-Rout X ^ ou t 

sin 2 i G 



8.8 x 10 10 sin~ 2 i M G 



(1) 



with -Rout = 7.5 kpc the outer radius of the rotating HI 
disc, v out = 225 km s" 1 the observed velocity of the disc 
at this distance, i the inclination of the disc (close to 90° 
for B2 0722+30) and G = 6.673 x 10" 11 m 3 kg -1 s -2 the 
gravitational constant. 



4.2 AGN classification 

A detailed continuum i mage of the rad io source in 
B2 0722+30 is presented bv lFanti et"afl (|l986l . also shown in 
Fig. 1 1 1 f) . From the continuum emission in our 21cm data we 
derive a total radio power of Pi.4ghz = 1.1 x 10 23 W Hz" 



for B2 0722+30, similar to the value found by iFanti et al.l 
{198$) . The o ptical V-band magnitude of B2 07 22+30 is 
m v = 1 5 04 iGonzalez-Serrano fe Carballdl2000l ). Follow- 
ing lOkd 0974), this corresponds to an optical flux den- 
sity of 3.5 x 10" 26 ergs s" 1 cm -2 Hz" 1 . The radio-to- 
opt ical flux ratio of B 2 072 2+30 is therefore -R1.4 ghz = 
41. iKellcrman n et al. ( 1989) defined R 5 ghz = 10 as the 
border between radio-quiet and radio-loud objects, which 
would correspond to roughly R ~ 19 at 1.4 GHz (see 
iKomossa etafl 120061 ). B2 0722+30 is therefore clearly a 
radio-loud object, although its total radio power lies in the 
transition region between that of Seyferts and radio galax- 
ies. Radio-loud Seyfert galaxies with total radio power and 
radio-to-optical flux ratio similar to that of B2 0722+30 do 
exist, but they are rare and their radio comp onent is gener- 
ally very compact (e.g. IKomossa et al. 2006, and references 
therein). An extr eme example is the nearby Seyfert 2 galaxy 
NGC 1068 (e.g. I Wilson fc Ulvestadl [19871 ) . which has a 1 
kpc-scale radio jet with a total radio power of the same or- 
der as B2 0722+30. In contrast to Seyfert galaxies, however, 
the radio source in B2 07 22+30 has a diameter of about 
13.6 kpc l|Fanti et al.lll987l ) and stretches outside the opti- 
cal body of the host galaxy. Although its linear size is still 
rather small compared to the majority of classical double- 
lobed radio sources, B2 0722+30 clearly contains a two-sided 
jet structur e (resembling so ca lled 'naked jets', as discussed 
in detail bv lParma et al.l 1987). The morphology of this two- 
sided jet extending over many kpc classifies B2 0722+30 as 
an FR-I radio source. 

In Sect. l3.3l we derived the emission-line ratios in the nu- 
clear region of B2 0722+30, which provide further insight in 
the properties of the AGN. The characteristic LINER classi- 
fication of the optical emission-line AGN is also in agreement 
with that of low-excitation FR-I radio galaxies. As already 
explained in Sect. 13.31 the emission- line ratios are incon- 
sistent with a typical Seyfert-type nucleus. Moreover, the 
[O III] luminosity is towards the low end of that found in 
Seyfert galaxies and the ratio of the 1.4 GHz radio power 
(Pi.4GHz = 1.1 x 10 23 W Hz" 1 ) over the [OKI] luminosity 
(£[0111] = 1-3 x 10 39 erg s _1 ) is about two orders of ma; 
nitude higher than the typical ratio for Seyferts IjWhittl 



ll985l ; lHo fc PengkOOll )! 3 ! Although intrinsic reddening likely 
underpredicts our derived [O III] luminosity in the nuclear 
region, we do not expect this effect to be significant enough 
to explain the relatively high radio-power/ [O III] ratio com- 
pared with that of Seyfert galaxies. Moreover, the emission 
lines in the diagnostic diagrams of IVeilleux fc Osterbrockl 
(1987) were specifically chosen to lie close together in wave- 
length in order to minimise reddening effects on the mea- 
sured line ratios. We therefore argue that intrinsic reddening 
effects do not alter our conclusions that B2 0722+30 has a 
LINER type optical AGN that is relatively weak compared 
to the radio source. 

From the combined radio continuum and optical 
emission-line properties of the AGN in B2 0722+30, we 
therefore classify this system as an FR-I radio galaxy rather 
than a genuine Seyfert galaxy. This implies that the low- 
excitation FR-I radio source would be hosted by the "wrong" 
kind of galaxy, given that classical radio sources are al- 
most always hosted by early-type galaxies rather than disc- 
dominated galaxies. 



4.3 The "wrong" kind of radio galaxy 

As discussed in detail in Sect. [T] classical double-lobed ra- 
dio sources are almost al ways hosted by early-type ga laxies, 
generally ellipticals (e.g. IVeron-Cettv fc VeronlhoOll ). The 
occurrence of an FR-I radio source in the disc-dominated 
galaxy B2 0722+30 is therefore exceptional, though not 
unique. Two other cases of disc-dominated radio galaxies 
have been studied in detail. The extended radio source 0313- 
192 is a classical FR-I source that is hosted b y a spiral 
galaxy (jLedlow et al.l 1998i . l200ll ; iKeel et al.l2006l ). NGC 612 
is an SO galaxy with a large star-forming HI disc, which 
hosts the powerful FR-I/FR-II radio source PK S 0131-36 
|Veron-Cettv fc Veron|[200ll ; lEmonts et al.ll200Sl ). Although 
the linear size of B2 0722+30 is significantly smaller than 
that of 0313-192 and NGC 612 and its total radio power is 
one (0313-192) to two (NGC 612) orders of magnitude lower, 
the fact that all three radio sources occur in a galaxy with 
a prominent disc and managed to escape the optical bound- 
aries of their hosts allows for a more detailed comparison 
between these systems. 

The fact that the occurrence of extended radio sources 
in disc-dominated host galaxies is extremely rare means 
that these three radio galaxies are excellent systems to in- 
vestigate in detail host galaxy properties or environmen- 
tal effects that may be important regarding the trigger- 
ing and/or evolu tion of radio-loud AGN. As explored by 
IKeel et al.l l|2006D and lEmonts et all l|2008l ). both 0313-192 
and NGC 612 share some striking similarities that could 
possibly be related to the occurrence of a classical double- 
lobed radio source in this "wrong" kind of host environment: 
NGC 612 and 0313-192 both contain a relatively luminous 
bulge compared with typical spiral galaxies, which generally 
reflec ts the presence of a relatively massive central black hole 



e.g. iKormendv fc Richstond Il995l ; iGebhardt et~ail 120001 ; 



iFerrarese fc Merrittll2000l) ; in both NGC 612 and 0313-192 



3 iHo fc Pend l l200ll) refer to P6cm, but based on references from 
the NASA/IPAC Extra-Galactic Database (NED), Pi. 4G Hz is of 
the same order of magnitude ELS -focm 

for B2 0722+30. 



© 2007 RAS, MNRAS 000.rTHT7l 



14 B. H. C. Emonts et al. 



the radio jet axis is oriented roughly perpendicular to the 
disc of the galaxy (at least within ~ 30°), which is likely the 
direction of least resistance from the ISM; in addition, both 
galaxies show indications that tidal encounters may have oc- 
curred, which could have resulted in enhanced fuelling of the 
central engine. 

For B2 0722+30 there is no clear evidence for the first 
two possibilities. B2 0722+30 appears to be a relatively 
small spiral galaxy (Sect. FO) , The total stellar population 
estimated from our SED modelling in the central kpc-scale 
region of the bulge - although uncertain - indicates that the 
bulge is not extremely massive. Although the power of the 
radio source in B2 0722+30 is one to two orders of magni- 
tude lower than that of 0313-192 and NGC 612 and there- 
fore not necessarily abnormal in the presence of a substan- 
tially smaller bulge region and lower mass black-hole, the 
average-mass bulge in B2 0722+30 does not provide an ob- 
vious explanation for why this particular disc-galaxy - com- 
pared with other spiral galaxies - contains an extended radio 
source. The radio jets in B2 0722+30 are also not aligned 
perpendicular to the disc of the host galaxy. The large-scale 
radio axis lies at a position angle (PA) of about 52° from 
the rotation axis of the disc. VLA-A re solution images of the 
radio continuum bv lFanti et al.l (| 1986h show a distorted jet- 
lobe structure, indicating that t he radio source changed d i- 
rection over its lifetime (see also Ide Ruiter fc Parmall 19841) . 
The jet direction in the inner kpc-scale region of the radio 
structure - tracing the more recently ejected jet material - 
lies even closer to the disc, at an angle of ~ 81° from the 
disc's rotation axis. 

It is interesting to note, however, that the orientation of 
the large-scale radio axis is aligned in the direction of the H I- 
rich MRK 1201/UGC 3841a pair. This is visualised in Fig. [IT] 
(left). Moreover, the orientation of the inner part of the radio 
jets is aligned in approximately the same direction as the H I 
bridge between B2 0722+30 and the MRK 1201/UGC 3841a 
pair. Although this could be merely a coincidence, we note 
that a similar alignment occurs in NGC 612, where the more 
recent radio-lobe material is aligned in the direction of an 
HI bridge towards a companion galaxy (Fig [11] - right). 
Perhaps a gas-rich interaction between B2 0722+30 and 
MRK 1201/UGC 3841a (as well as between NGC 612 and 
its companion) somehow cleared a path for the radio source 
to escape the host galaxy with minor resistance from the 
ISM, possibly in combination with an induced precession of 
the accretion-disc/jet-axis in this direction. In a future pa- 
per, we will investigate in more detail whether an alignment 
between the radio axis and the direction of Hl-rich com- 
panions or tidal material occurs also in other nearby radio 
galaxies. 

4.4 Interactions and radio source triggering 

Interactions likely occurred between B2 0722+30 and nearby 
companions, which could be related to the triggering of the 
radio source in B2 0722+30 and/or to the fact that the ra- 
dio source managed to escape the optical boundaries of the 
host galaxy. The interacting pair MRK 1201/UGC 3841a 
west of B2 0722+30 comprises galaxies that appear to 
be in ongoing interaction, based on their HI properties. 
The faint HI bridge in between B2 0722+30 and the 
MRK 1201/UGC 3841a system suggests that B2 0722+30 is 



involved in the interaction. The distribution of H I gas in the 
MRK 1201/UGC 3841a pair in the vicinity of B2 0722+30 
is very similar to that of the NGC 5058 group, where an 
Hl-bridge bet ween two galaxies winds ar ound one side of a 
shell elliptical dGuhathakurta et a.1.11 1990h . Further tentative 
tidal features that may indicate that B2 0722+30 has been 
involved in recent interactions are the HI extension at the 
NW-tip of the disc of B2 0722+30 in the direction of a pos- 
sible faint companion ("XI" in Fig. [1} and - if confirmed - 
the fainter HI tail further north. 

It is likely that galaxy interactions perturbed the gas in 
the central part of the disc of B2 0722+30. This could have 
resulted in a loss of angular momentum of the gas so that 
it may have been deposited on to t he central black- hole and 
triggered t he radio so u rce (s ee e.g. Lin et all 19881). Recent 
studies by iKuo et all l|2008h and iTang et al.l (|2008l ) show 
that Hl-rich galaxy interactions are generally also associ- 
ated with Seyfert galaxies, much more than with their non- 
active counterparts. We would like to note, however, that 
in an ongoing HI study of a larger sample of nearby low- 
power radio sources (all hosted by early-type galaxies) we 
do not see evidence in HI for ongoing gas-rich interactions 
associated with extende d FR-I radio source s (Emonts et a l. 
in preparation; see also lEmonts et aT]|2007l ; lEmontsll2006l ). 
In fact, recent studies suggest that the triggering mecha- 
nism for most FR-I radio sources - or low-excitation radio- 
AGN in general - may be related to Bondi accretion of hot 
circum-galactic g a s on to the central bl a ck hole (see e.g. 
Allen et alj 120061 ; lHardcastle et all 120071 ; iBaldi fc Capettil 
20081 ). Although our data cannot rule out a similar accre- 
tion mechanism for B2 0722+30, the results presented in 
this paper suggest that - at least in certain systems and un- 
der the right conditions - low-power FR-I radio-AGN can 
occur in systems that are undergoing gas-rich galaxy in- 
teractions. Whether or not the interactions associated with 
B2 0722+30 are directly responsible for the triggering of 
the AGN through fuelling of the black hole with low angu- 
lar momentum gas needs to be investigated by studying the 
kinematics of the gas in the central region at higher resolu- 
tion. Further investigation into this matter is important for 
understanding the triggering mechanism of various types of 
AGN. 



5 CONCLUSIONS 

B2 0722+30 is a disc galaxy that hosts a classical double- 
lobed radio source. Using HI radio synthesis observations, 
deep optical imaging, stellar population synthesis modelling 
and emission-line diagnostics we investigated in detail the 
properties of the host galaxy and radio-loud AGN as well 
as the environment of B2 0722+30. The most important 
conclusions from this paper are: 

• B2 0722+30 contains a regular, rotating disc of both 
neutral hydrogen and ionized emission- line gas. A reddened 
young stellar population is present throughout this disc and 
emission- line ratios are characteristic of HII regions, which 
indicates that star formation is ongoing in the disc. The 
high FIR luminosity of B2 0722+30 may originate from 
dust heating by the young stellar population; 
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B2 0722+30 NGC612 




Figure 11. Left: Apparent alignment of the radio jets/lobes of B2 0722+30 in the direction of the interacting galaxy pair 
MRK 1201/UGC 3841 a. The contours i n the main panel show the total intensity HI map, while those in the inset show the radio 
continuum (taken from lFanti et al]|l98rj) . The solid line shows the original direction of the radio lobes (PA = 97°), which is towards 
the extended optical tail of the distorted and elongated disc galaxy UGC 3841a. In the inner region (at kpc-scale distance from the 
nucleus), more recently ejected radio plasma apparently propagates roughly along PA = 124°, which is in the direction of the HI bridge 
in between B2 0722+30 and the MRK 1201/UGC 3841a pair (dashed line). Right: a similar alignment is seen for NGC 612, where the 
most re cently ejected radio plasma propagates in the direction of a large-scale HI bridge that extends towards a gas-rich companion 
galaxy ( Emon ts et al 1 l2008t). Agai n, the cont ours in the main panel show the total intensity HI map, those in the inset show the radio 
continuum (taken from iMorganti et al.|[l993h and the solid line indicates the position-angle of the radio lobe. 



• The AGN characteristics classify B2 0722+30 as an 
FR-I radio galaxy rather than a Seyfert galaxy; the radio 
source has an FR-I morphology, the optical AGN shows 
emission-line ratios characteristic of a LINER nucleus and 
the [O III] luminosity, in particular when compared to the 
radio power, is relatively low; 

• The presence of a classical double-lobed radio source 
in a disc-dominated host galaxy is extremely rare in the 
nearby Universe. Radio galaxies like B2 0722+30 therefore 
provide an excellent opportunity to identify host galaxy 
properties and environmental effects that could be related 
to the triggering mechanism of radio- loud AGN; 

• B2 0722+30 has an Hi-rich environment, with several 
gas-rich galaxies that appear to be interacting. We argue 
that interactions between B2 0722+30 and nearby systems 
could be related to the triggering of the radio source. How- 
ever, more detailed kinematical studies of the gas in the nu- 
clear region of B2 0722+30 are necessary to confirm this and 
to better explain the physical processes involved. The align- 
ment of the radio axis in the direction of the gas-rich galaxy 
pair MRK 1201/UGC 3841a offers the interesting possibil- 
ity that an interaction with this pair may be related to the 
fact that the radio source managed to escape the optical 
boundaries of the host galaxy. 

B2 0722+30 is part of a complete sample of nearby B2 
radio galax ies for which we studied the large-s cale H I prop- 
erties (see lEmontsI [200rj ; lEmonts etaLI 120071 ). A detailed 
analysis of the complete sample will be given in a forth- 
coming paper. 
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